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Abstract 


We simulate the formation of cosmic strings at the zeros of a complex 
Gaussian field with a power spectrum P{k) oc kP, specifically addressing the 
issue of the fraction of length in infinite strings. We make two improvements 
over previous simulations: we include a non-zero random background field in 
our box to simulate the effect of long-wavelength modes, and we examine the 
effects of smoothing the field on small scales. The inclusion of the background 
field significantly reduces the fraction of length in infinite strings for n < —2. 
Our results are consistent with the possibility that infinite strings disappear 
at some n = Uc in the range —3 < nc < —2.2, although we cannot rule out 
Uc = —3, in which case infinite strings would disappear only at the point 
where the mean string density goes to zero. We present an analytic argument 
which suggests the latter case. Smoothing on small scales eliminates closed 
loops on the order of the lattice cell size and leads to a “lattice-free” estimate 
of the infinite string fraction. As expected, this fraction depends on the type 
of window function used for smoothing. 


Typeset using REVTgX 


1 


I. INTRODUCTION 


Cosmic strings are effectively one-dimensional topological defects which may form at a 
phase transition in the early universe (see Ref. [Q for a review). Although much of the early 
interest in cosmic strings has centered on the possibility that they might have served as seeds 
for the formation of large-scale structure, cosmic strings are interesting physical objects in 
any case, and they have analogues in the study of condensed matter [Q. 

Any investigation of the evolution and cosmological consequences of cosmic strings must 
begin with the study of the initial cosmic string conhguration, a study which was hrst 
undertaken by Vachaspati and Vilenkin |P], and re-examined by many others ( |@] - 1^). 
Although subsequent cosmic string evolution will erase many of the details of the initial 
conhguration, one fundamental property of the initial conditions is crucial to the subsequent 
evolution: the existence of inhnite strings. Without the existence of inhnite strings, the 
cosmic strings produced at the phase transition may all decay via gravitational radiation 
long before they can have any interesting cosmological effects. Vachaspati and Vilenkin 
found in their simulation that roughly 80% of the string was in the form of inhnite strings 
0]. Using a very diherent type of simulation for the string formation process, Borrill P] 
claimed that this fraction was, in fact, zero, while Robinson and Yates [Q, in a study of the 
dependence of this fraction on the power spectrum of the initial held, argued that for power 
spectra of the form P{k) oc fc”, the inhnite string fraction f^o drops to zero for n < —2. In 
fact, this result is not obvious from their simulations; it is based on htting f^o to an analytic 
function for n > —2. 

In this paper, we extend the simulations of Robinson and Yates in two ways. First, 
we include the ehects of long-wavelength modes which are absent from earlier simulations. 
Second, we examine the ehect of smoothing the initial held to remove lattice ehects. We 
hnd that the former change has a dramatic ehect on foo, sharply reducing f^o for n < —2. 
Smoothing also ahects /oo, leading to “lattice-free” estimates of the fraction in inhnite 
strings. In the next section we present our numerical results, and in Section 3 we discuss 
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briefly our main conclusions. 


II. NUMERICAL SIMULATIONS 


With one exception [§], simulations of cosmic string formation typically make use of a 
lattice. A value of the field 0 associated with the string is assigned to each of the cells of 
the lattice, and the location of the cosmic string is then identified with edges of the lattice 
around which the field winds through 360". The original simulations of Vachaspati and 
Vilenkin [Q were performed on a cubic lattice with no correlations between values of the 
field in different cells. Subsequent researchers investigated changing the probability of string 
formation by “biasing” the distribution of the values of the field p, P], 0, by allowing the 
field to be divided into domains of variable size @], or by allowing long-range correlations 
between the field values [Q. All simulations on a cubic lattice suffer from the problem that 
an ambiguity exists with regard to the string assignments at the vertices of the cubes; this 
problem can be eliminated by going to more exotic lattices, such as the tetrakaidekahedral 
lattice, for which only four edges meet at every vertex i, 0- 

In this paper, we simulate the formation of cosmic strings using a cubic lattice and a 
complex Gaussian field 0 with long-range correlations, where the strings are taken to he 
along the zeros of 0. A Gaussian field is completely characterized by its power spectrum, 
defined by 

P{k) = I d\ e**^'’'{0(x)0(x + r)). (1) 

We take 0 to have a power-law power spectrum: 


P{k) oc iV 


( 2 ) 


Models of this sort were first examined by Vishniac, Olive, and Seckel |^, who proposed 
that cosmic strings with n = —3 could be produced due to quantum fluctuations of the field 
0 during inflation. Vishniac et ah showed that the mean string length per unit volume, 
L/V is given by JTO 
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( 3 ) 


where 

_ JP{k)k^dk 
^ ~ JP{k)k^dk' 

Note that some sort of smoothing or cutoff is required for convergence at large k, but this is 
provided automatically by the lattice cut-off in numerical simulations. As n —> —3, equation 
(§) gives LjV 0; this result is conhrmed by our numerical simulations. Physically, what 
happens is that within arbitrarily large regions of space the held 0 has Re(0) > 0 (for 
example). An expression similar to equation (^) is also given in reference 

In our model, we assign values of cf) to the sites on a cubic periodic lattice, and the string 
is assumed to lie along the location of the zeros of this complex held. This model, as we 
have described it, is identical to that of reference [^, but we make two improvements in the 
model, both of which have to do with the limitations in dynamic range inherent in such a 
simulation. 

Consider hrst the largest scales. Because our simulation volume is hnite, the simulation 
loses power on all scales larger than the box size, a problem which was noted by Robinson 
and Yates Q. In the simulations, this corresponds to the fact that the mean value of 0 
averaged over the entire box is zero. This problem can be resolved by adding a random 
uniform background held 0b to every cell in the box; this held represents the contribution to 
the held value from all Fourier modes which are larger than the box size. The variance of 0b 
can be determined from equation (|^), and in our simulations it is chosen to have a Gaussian 
distribution with this variance. 

Our second addition to the simulation involves the behavior of the held on small scales. 
The use of a lattice to simulate string formation is obviously unphysical, and numerical 
simulations clearly show that foo depends on the lattice being used. (Compare, e.g., the 
results of reference for a cubic lattice with those of [|] and 0 for a tetrakaidekahedral 
lattice.) On the other hand, we do expect the held 0 to be correlated on the smallest scales. 


{k^ 
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leading to some sort of domain strnctnre. We resolve this problem by smoothing the held 
on a scale larger than the cell size to eliminate lattice effects. 

Our simulations were performed on a cubic lattice of size 128^. A Gaussian, complex¬ 
valued random held </> having the power-law power spectrum given in equation (|^) was set 
down on the lattice. The strings were identihed as vortices of the held 0. In tracing a 
closed path around the four cells bounding each edge of the lattice, the held was assumed to 
change values by moving in the shorter of the two possible directions in the complex plane. 
A cosmic string was then placed along an edge if the held traced out a 360° circle in the 
complex plane as a path was traced around the four cells bounding that edge. We dehne an 
inhnite string as one which crosses all the way from one end of the box to the other in at 
least one of the three directions. Note that this dehnition dihers from that used by Robinson 
and Yates 0, who used a cutoh in the string length. 

Consider hrst the simplest case, for which we have no smoothing, and the mean held 
value is set to zero. This corresponds exactly to the simulations of reference [^j. To hnd the 
value of /oo and its variance, we performed 32 simulations for each value of n, grouped into 
eight groups of four simulations. Within each group of four, we derived an average value for 
/oo by dividing the combined length of inhnite strings in all four simulations by the total 
length of string in that group. This procedure corresponds to considering each of the four 
simulations in a group as sampling a diherent region of space. We then averaged /oo for the 
eight groups of simulations to derive a hnal mean and variance. [Note that in deriving 
Figures 3-6, we used 4 groups of 4 simulations rather than 8]. 

In Figure la, we present our results for no smoothing with zero mean held. (All error 
bars are l-a). These results agree closely with those of reference 0, as they should. This 
indicates that /oo is relatively insensitive to the exact dehnition of inhnite string, since this 
is the only diherence between our simulations and those of reference [Q. 

Now we repeat the simulations, but add a random uniform background held 0b to the 
entire box. This background held has a Gaussian distribution with a variance which can 
be determined, in principle, from equation (|^). In practice, we use a simpler method: we 
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embed our simulation volume in a larger 128^ box with the same power spectrum, in which 
our entire simulation volume occupies a single cell of the larger box. The mean held (j)b in 
the simulation volume is then simply the value of 0 in a single (randomly-chosen) cell of 
the larger box. Note that this method still causes a loss of power at the largest scales (i.e., 
Fourier modes with wavelengths longer than the size of the larger box), but it provides an 
effective dynamic range of 128^ = 16, 384, which represents a considerable improvement over 
the zero mean held case. 

In Figure lb, we show our results for this case. The greatest diherence between Figures 
la and lb occurs for n < —2: the addition of the background held produces a larger 
variance between runs, and it reduces foo sharply. The graph in Figure lb is consistent 
with the suggestion of reference 0 that goes to zero when n is less than some critical 
value Uc, where > —3; i.e., the average string density is still non-zero when the inhnite 
strings disappear. In reference [^], it was claimed that ric = —2, but this is inconsistent 
with our results. If we exclude points for which f^o is more than 3-a from zero, then we hnd 
He < —2.2. Although our results are consistent with the possibility that ric > —3, they do 
not rule out the possibility that ric = ~3 (i.e., the inhnite strings do not disappear until all 
the strings disappear). 

In Figure 2, we plot the total length in closed loops and inhnite strings as a function of n 
with the inclusion of the background held. (All of the results presented in this paper, with 
the exception of Figure la, include this background held). Although the length in inhnite 
strings is signihcantly reduced below that obtained without a background held (e.g., see 0), 
we cannot tell with certainty at which value of n the inhnite strings disappear. 

In Figure 3, we show the size distribution of closed loops for two representative cases 
(n = 0 and n = —2.8) when the background held is included. In both cases, N{1) follows a 
power law in /, but the slope changes noticeably with n. Fitting the points with 10 < / < 100, 
we hnd that 


N{1) OC /- 2 - 62 ± 0 . 01 ^ 


( 5 ) 
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N{1) OC /-2-35±0.04^ _ _2_g)_ 

These fits are shown as dashed lines in Figure 3. Neither of these is consistent with the 
frequently assumed behavior N{1) oc but even the original {n = 0) simulations 

of Vachaspati and Vilenkin [Q found that N{1) oc /- 2 - 6 ±o.i_ 

Now consider the effect of smoothing the (j) held on small scales. For a window function 
W{r), the smoothed held (/>s(x) is the convolution of 0(x) with hF(r); 

0s(x) = y 0(x + r)hF(r). (7) 

The ehect of smoothing is to reduce the magnitude of the small-scale huctuations by aver¬ 
aging them out over the window function. We calculate foo when 0 is smoothed with three 
diherent window functions, the spherical Gaussian window function, 

hF(r) = exp(—r^/2rQ), (8) 


the spherical tophat. 


fF(r) = 1 (r < To), 

W(r) = 0 (r > To), (9) 


and a sharp cut-oh in fc-space, which corresponds, in physical space, to the smoothing 


sin(r/ro) cos(r/ro) 

(r/ro)3 (r/ro)2 


( 10 ) 


where we then normalize each window function to give fW(r)d^r = 1. 

We now examine the variation of foo with rg. We consider only the case n = 0 because 
it resembles the most likely scenario for string formation (i.e., no long-range correlations in 
the 0 held), and for the range of n values we have examined, it has the most short-range 
power and should therefore show the largest sensitivity to smoothing. Our results are shown 
in Figure 4a for the spherical Gaussian, in Figure 4b for the spherical tophat, and in Figure 
4c for the sharp cut-off in fc-space. 
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Note that ro represents something different for each of our three window functions, so it is 
meaningless to compare foo for the same value of tq in each figure. However, the important 
point is whether there exists a range of values of ro for a given window function which 
produces a relatively constant set of values for /oo- For the case of Gaussian smoothing, 
foo decreases with tq but eventually “plateaus” for tq > 1.5, suggesting that the value of 
foo in this region is the “true” value of foo for Gaussian smoothing. Averaging foo for 
1-5 < ro < 2.5, we obtain foo = 0.71 ± 0.01. As ro is increased, the size of the smoothing 
volume grows compared to the size of the box, and we see larger fluctuations from one run 
to the next. However, the rms fluctuations in foo remain within 10% of our mean value 
(/oo = 0.71) over the entire range of ro values we have examined (ro < 10). An average of 
foo over the range 1.5 < ro < 10 gives foo = 0.72 ± 0.04. 

For spherical tophat smoothing, we again see that foo decreases with ro, but no obvious 
plateau is visible in Figure 2b. In fact, tophat smoothing fails to eliminate lattice effects. 
If we substitute the smoothed version of the power spectrum into equation (|[), we find 
that for the spherical tophat window function with n = 0, the integral in the numerator in 
equation (^) fails to converge at fc —> cx), indicating that it is the cubic lattice, rather than 
the spherical tophat, which provides the cutoff in this case. This is confirmed by the loop 
distribution for the spherical cutoff; unlike our other two window functions, the spherical 
tophat produces a loop distribution which retains a large number of loops with size of order 
the cell size. 

The sharp fc-space cut-off shows qualitatively similar behavior to the Gaussian window 
function. The value of foo decreases slightly as a function of tq, and reaches a plateau for 
ro > 1. Averaging foo for 1 < ro < 3 and 1 < ro < 10, we get foo = 0.82 ± 0.01 and 
foo = 0.85 ± 0.05, respectively. 

In Figure 5, we demonstrate the effect of smoothing on the size distribution of the loops, 
by showing N{1) for the closed loops with Gaussian smoothing and rg = 3. Note that 
smoothing does not eliminate the closed loops smaller than the size of the window function; 
rather, it reduces the number of loops to a constant as a function of I for small 1. Of course. 
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the total length in closed loops is then dominated by loops with size near the point at which 
the power law behavior begins (in this case, I ~ 30). 


Smoothing allows us to test the analytic formula for the total string density given in 
equations and (^. For the power spectrum given by equation (^) smoothed with the 
Gaussian window function given by equation (^, the expression for LjV reduces to 



( 11 ) 


(Note that in this case, it is obvious that LjV —>■ 0 as n —*• —3). To compare with our 


simulations on the lattice, LjV must be multiplied by the number of cells (128^ in our case) 


times the number of independent lattice edges per cell (3 for a cubic lattice). The expected 
string length derived in this way from equation (|TT|) is compared to our simulation results 
in Figure 6, for a smoothing length tq = 3. The agreement is excellent, within the expected 
statistical fluctuations. Whether one views this as a conhrmation of the analytic expression 
or of the validity of our numerical simulations is a matter of personal preference. 


III. DISCUSSION 


Our results are consistent with the possibility that a transition occurs in the string 
network at some value n = Uc, with —3 < < —2.2; when n < ric the inhnite strings 


disappear from the network, leaving only closed loops. Robinson and Yates 0 argued that 


ric = —2, but this value is not consistent with our results. Furthermore, we cannot rule out 
the possibility that Uc = —3, i.e., the inhnite strings do not disappear until the mean string 
density goes to zero. 

In fact, the following argument suggests that inhnite strings should not disappear for 
n > —3. The real and imaginary parts of 0 are independent Gaussian helds, and the zeros 
of each of these helds form a set of two-dimensional surfaces. Gonsider hrst the surfaces 
dehned by Re(0) = 0. Since the volumes of space occupied by the regions with positive and 
negative Re(0) are equal, we expect both regions to percolate to an inhnite distance. Hence, 
the boundaries dividing these regions should contain at least some inhnite surfaces. This 
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argument holds for both the surfaces dehned by Re(0) = 0 and lm(0) = 0. The intersection 
of these two sets of surfaces gives us the location of the cosmic strings. Since inhnite surfaces 
must exist in both sets of surfaces, it appears that there must also be inhnite cosmic strings. 
Note that it is possible to imagine rather arcane distributions of the helds which violate 
this argument. For example, the regions with Re(0) > 0 and Re(0) < 0 could be nested 
inside of each other in larger and large hnite volumes, producing a fractal distribution with 
arbitrarily large but hnite surfaces of Re(0) = 0. However, it seems unlikely that a Gaussian 
held could lead to such a distribution. The crucial point in this argument is the fact that the 
distribution is symmetric with respect to positive and negative values of Re(0) and lm(0); 
if we relax this assumption and “bias” the distribution, the argument no longer holds. In 
fact, inhnite strings are observed to disappear in simulations with such a “bias” i, i, 0- 

The aim of our simulations with a smoothed held 0 was to obtain a lattice-independent 
estimate of foo- We considered only the white noise spectrum, n = 0, and found foo ~ 0.7 
for spherical Gaussian smoothing and foo ~ 0.8 for a sharp cut-oh in /c-space. These values 
are comparable to those obtained in earlier lattice simulations without smoothing The 
variation of foo for diherent choices of smoothing is not surprising, since, for example, the 
total string density L/V in equations and (^) clearly depends on the short-wavelength 
behavior of the spectrum P{k) (which is ahected by the smoothing). This variation is again 
comparable to the variation between the values of foo in simulations without smoothing on 
different types of lattice @], 0, 0]. 

We hnd no evidence supporting the hypothesis 0 that the presence of inhnite strings is 
due entirely to the lattice and that they should disappear in lattice-free simulations. The 
smoothing length we used was sufficiently large to ensure that the loops making the largest 
contribution to the total string length had sizes much greater than the lattice cutoff, so that 
the lattice effects were minimal. Still, we found a substantial fraction of the total length to 
be in inhnite strings. 

We are currently investigating the formation of domain walls and monopoles with cor¬ 
related helds. Given that a parallel literature on this subject exists in condensed matter 
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physics [^], these results may have applications beyond the purely cosmological. It is con¬ 
ceivable, for example, that the decline of the inhnite string density with the decrease of the 
spectral index n can be tested experimentally. What one needs is a condensed matter system 
(such as liquid He^) in which linear defects are formed at a second-order phase transition. 
Defect formation can then be observed by a rapid temperature (or pressure) quench from 
above to below the transition point |ll^,[l^. Near the critical temperature Tc, the order 
parameter develops long-range fluctuations. At T = Tc, the fluctuation spectrum is a power 
law of the form (|^) with n = —2 -1- rj, where the critical exponent r] is typically a small 
number, r] ^ 0.05 [jl^ {rj ^ 0.05 for He^). If the system is allowed to equilibrate very close 
to the critical point and is then rapidly quenched to subcritical temperatures, one can expect 
the length in inhnite strings to be suppressed compared to a quench from a temperature 
well above Tc (where the huctuation spectrum is close to n = 0). Our Fig. 2 suggests a 
suppression roughly by an order of magnitude, while the value of Uc = —2 conjectured by 
Robertson and Yates 0 would give a much more dramatic suppression. It should be noted 
that a realistic quench is a rather complicated process, and its outcome can depend on a 
variety of physical effects (for a recent discussion see P^JT5[]). 
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FIGURE CAPTIONS 


Figure 1: The fraction of total string length in the form of inhnite strings, /oo, as a function 
of n, where P{k) oc /c”, and P{k) is the power spectrum of the complex Gaussian held 
which gives rise to the strings, where (a) the held in the box is forced to have zero mean, 
and (b) a background non-zero mean held is added to the box to simulate the ehects of 
long-wavelength modes. 

Figure 2: The total length L in closed loops (open squares) and inhnite strings (crosses) as 
a function of n, where P{k) oc and P{k) is the power spectrum of the complex Gaussian 
held which gives rise to the strings. A background non-zero mean held is added to the box 
to simulate the ehect of long-wavelength modes. 

Figure 3: The number of loops N with a given length / for (a) n = 0 and (b) n = —2.8, 
where P{k) oc /c”, and P{k) is the power spectrum of the complex Gaussian held which gives 
rise to the strings. A background non-zero mean held is added to the box to simulate the 
ehect of long-wavelength modes. The dashed line is the best-ht power law in each case. 

Figure 4: The fraction of total string length in the form of inhnite strings, /oo, as a function 
of smoothing length tq for (a) Gaussian smoothing, (b) spherical tophat smoothing, (c) 
a sharp /c-space cutoh, where the power law index is n = 0 (no held correlations), and a 
background non-zero mean held is added to the box to simulate the ehect of long-wavelength 
modes. 

Figure 5: The number of loops N with a given length I for a held with n = 0 (no held 
correlations) smoothed with a Gaussian window function with smoothing length tq = 3. A 
background non-zero mean held is added to the box to simulate the ehect of long-wavelength 
modes. 
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Figure 6: The total string length L as a function of n, where P{k) oc /c”, and P{k) is 
the power spectrum of the complex Gaussian held which gives rise to the strings. The 
points with error bars are the results of our simulation, and the solid line is the analytic 
prediction. A background non-zero mean held is added to the box to simulate the ehect of 
long-wavelength modes. 
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